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We define a complete set of scalar master integrals (Mis) for massive 2-loop QED Bhabha scattering. Among 
others, there are thirty three 2-loop box type Mis. Five of them have been published in (semi-)analytical form, 
one is determined here, the rest remains to be calculated. Further, the last four so far unknown 2-loop 3-point 
Mis are identified and also computed here. 



1. INTRODUCTION 

The determination of luminosity at high energy 
colliders is a backbone linking theoretical calcu- 
lations with experimental data. Bhabha scatter- 
ing, e~^e~ — > e"^e~(ri7), is the process that al- 
lows to accomplish this task. So far the complete 
NNLO corrections in massive QED have not been 
computed. There are several reasons to fill this 
gap. First, the future experimental accuracy is 
expected to be of the order of 10"'* |T] . Obviously, 
the theoretical error should be smaller. Second, 
the current Monte Carlo generators for real pho- 
ton radiation assume massive external fermions. 
Analytical results of NNLO corrections in mas- 
sive QED are clearly needed there for consistency. 
There are also technical reasons to perform the 
calculation. First, as argued in pi, 2-loop QED 
results are a useful testing ground for 2-loop QCD 
physics containing more than one kinematic in- 
variant. Second, the full analytic results can be 
easily approximated to various kinematic regions; 
this is not possible for massless results. Finally, 
the complete calculation of virtual 2-loop correc- 
tions in massive QED is a prelude to the deter- 
mination of analogous effects within the Standard 
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Figure 1. The 2-loop self energy diagrams. 



In this contribution, we report on first results 
of our calculation of virtual 2-loop corrections in 
massive QED. We have determined a complete set 
of scalar master integrals to be calculated for the 
problem at hand, including 33 2-loop box master 
integrals^. Subsequently, we have calculated all 2- 
and 3-point Mis, four of which were still lacking 
in the literature. We also started the calculation 
of boxes. 



^More detailed information, including the Mis computed 
so far, can be found at the webpage bhabha.html in |5]. 



2. 2 LOOP DIAGRAMS AND IDENTI- 
FICATION OF Mis 

In Figs, n to 131 we show the set of 2-loop self 
energy, vertex and box diagrams for which the 
scalar Mis are needed. 







Figure 2. The 2-loop vertex diagrams. 
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Figure 3. The 2-loop box diagrams. 

Following the Laporta-Remiddi approach |6I7| , 
we have determined a set of scalar Mis for the vir- 
tual 2-loop corrections to Bhabha scattering. We 
use the C-I--I- library DiaGen/IdSolver'^, which 
determines and solves an appropriate set of alge- 
braic equations, derived with integration by parts 
(IBP) ^ni and Lorentz invariance (LI) identities 
|11| . The latter have been useful for algorithmic 
optimization but did not reduce the number of 



^DiaGen has already been used in several other projects 
I8I9I . IdSolver is a new software package by M.C. We are 
also using Fermat ^, FORM |31, Maple, Mathematica. 



Mis. In Tabled we give a list of the net numbers 
of 2-loop master integrals needed for the evalu- 
ation of the various 2-loop vertex and box dia- 
grams. Obviously, the 1-loop masters add up. A 
complete list of 2-loop Mis is given in Figs. 01 to 

m 
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Figure 4. The six 2-loop 2-point Mis. The multi- 
plication factor stands for the number of masters 
for the given diagram when different from 1. 

Strictly speaking, we show prototypes, i.e. 
topologies with masses assigned to the lines. 
Some of the prototype diagrams are multiplied 
by an integer. This indicates that masters ap- 
pear with dotted lines (higher powers of propaga- 
tors)^. As an example, we give the master list for 
the prototype V512m (i.e. V-ertex with 5/-ines, 
among them ^m-assive). All the other masters 
may be found in 0. There, we also give tables 
indicating which Mis contribute to SE1-SE4, VI- 
V5, B1-B6. 

We have started to solve the Mis in a system- 
atic way with the differential equation approach 
proposed in 1^ and developed in T^^. With 
(pi + P2) = {P3 + Pa) and pf = ml, we use 
s = {pi+p2y, t = (pi—ps)^- The analytical 
results depend on the variables x and y, 



V-s + 4-h 



(1) 



and y is obtained by replacing s by t. We set 
everywhere the electron mass to unity, rrie ~ 1. 

*One might also use irreducible numerators. However, this 
introduces an explicit dependence of the master on the 
specific momentum distribution chosen. 
^We put the results obtained so far in the Mathematica file 
MastersBhabha.m in |^. Recalculating some of the known 
results, we found some misprints in Eqs. (80), (85), (121), 
(128), (131), (144) in [H and in Eqs. (36), (38), (39), (B. 9) 
in 1151 for the corresponding masters. We have been in- 
formed by the authors that an Erratum is in preparation. 



Table 1 

The number of 2-loop master integrals needed to calculate the 2-loop vertex diagrams (Fig. |2) and box 

diagrams (Fig. ^ 



Diagram 


VI 


V2 


V3 


V4 


V5 


Bl 


B2 


B3 


B4 


B5 


B6 


2-point Mis 
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3-point Mis 
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10 
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11 


13 


10 
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Box type Mis 
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Figure 5. The nineteen 2-loop 3-point Mis (the 
multiplication factors as in Fig. ^) 
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Figure 6. The thirty three 2-loop box Mis (the 
multiplication factors as in Fig. 0|) 



Our normalization of the momentum integrals 
is chosen such that the one-loop tadpole is {d = 
4-2e): 

d'^q 



Tlllm = 



(TreT^)"^ 



1 



92-1 

4 



(2) 



3. NEW 2-LOOP 3-POINT MASTERS 

Most of the nineteen 2-loop 3-point Mis needed 
for the massive Bhabha process have been cal- 
culated already: V412m (2 masters), V413m (2 
masters), V414m (2 masters), V513m (2 masters), 
V614ml (2 masters) in ^.■, V411m (4 masters) and 
V614iii2 in ^B]. The last four masters completing 



the list, needed for box diagrams Bl and B3, are 
of the type V512m (see Fig. [Tjl. We obtain the 
following analytical expressions: 



V512ml 



-^a + C2H{0,0,x) 



(-1 + .T2) 

-2C2i?(l,0,x) + 2iJ(0,0,0,l,a;) 

+H{0, 0, 1, 0, x) + 2iJ(0, 0, 1, 1, x) 

+2H{0, 1, 0, 0, x) - 3iJ(l, 0, 0, 0, x) 

+2H{1, 0, 0, 1, x) + 2iJ(l, 0, 1, 0, x) 

-AH{1, 0, 1, 1, x) - 4C3i?(0, x) - 2C3i?(l, x) 

+0{e), (3) 



V512mld = 



2x 



(-1 + .x)2(1 + .t) _' 



xH{Q,0,x) 



+ (-l + a;)[2C2 + ff(0,l,2:)] 



2x 



+ 0{e), (4) 



V512m2 = . , ,, 8C4 + 2C2i?(0, 0, x) 
(-1 + x^) [ 

-4C2i?(0, 1, x) + H{Q, 0, 0, 0, x) 



+2i?(0, 0,0,1,2;) 



OW, 



(5) 



V512m2d 



£2 4(-l + a:)2 



e 2(-l + a;)2 



-l + iJ(0,x)+2i7(l,j;) 



-4 - I5C2 - 4a; 



4(-l + x)2(l + j;) 

-C22; + 4(1 + a;) [H{Q, x) + 2i/(l, x) 

-2(a; + 3) [iJ(0, 0, x) + 2iJ(0, 1, x) 



-8(1 + x) [i/(l, 0, x) + 2H{1, 1, x)] 



0(6). (6) 



The ^(0, . . . ,x) etc. are harmonic polyloga- 
rithnis as introduced in |17| . We have performed 
a variety of cross-checks on the results, e.g. nu- 
merical evaluations at fixed kinematical points in 
the Euclidean region using sector decomposition 
following the approach described in "18J . The nu- 
merical integrations are done with CUBA J^. See 
also 1201 for further cross checks. 

4. A NEW MASSIVE 2 LOOP BOX MI 

As mentioned, there are thirty three 2-loop box 
Mis needed. By now, few of them are known an- 
alytically. For diagram Bl, the masters B714ml 
and B714mlN are given in |2J|; for diagram B2, 
the master B714m2 is known as two-dimensional 
integral representation [22]; for diagram B3, the 
leading divergent part of master B714m3 is pub- 
lished [22]. For diagram B5, the two double box 
masters, B514m and B514md, have been derived re- 
cently with the restriction to Nf = 1 (one flavor) 
in 15 5 ^ 

A further double box master integral, B512ml 
(see Fig. E)), is given here; it contributes to dia- 
grams Bl and B3 and is the solution of the fol- 
lowing differential equation: 
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Figure 7. Four Mis of type V512m. 



S-— B512ml(s,t) = 

OS 

\ + s^ -2t + s(-6 + t + et)' 
(-4 + s)(-4 + s + t) 

"2-9e + 9e2 



€{-A + s)t 
(-l + 3e)(-4 + t) 



(-4 + s)(-4 + s + t)_ 
s(l-3e) 



(-4 + s)(-4 + s + f) 



B512ml(s,t) 
SE310m(t) 
V411m2(t) 
V412ml(s). (7) 



To get this equation we use a simple represen- 
tation for differential operators. In the s-channel 
(for the t-channel operator: p2 *-^ ~Ps,) 



.^ = i|(p5'+P^)-^^^'^^^ 



ds 2 1^'^' ' ""^^ ' s + t-Am?\ dpt^' 

The boundary condition^ is fixed at a; = 1. 
The solution is: 



(8) 



1 
B512ml = — 



(-l + a;2) 



H{0,x) 



e (-1-Ha;2) 



C2-H{Q,x)[2 + H{0,y) 



®We have also used the point x = — 1 or on-shell conditions 
with cross checks with the 0NSHELL2 package 1231 . 



-2H{1, y)] + 2H(-1, 0, x) - H{0, 0, x) 



-G{^y,x)[3C2 



(-1 + :e2) 

+H{0, 0, y) + 27f (0, 1, y)] + G{-l/y, x) [6(2 

+F(0,0,j/) + 2i/(0,l,y)]-2C2 

+2G(-l/y, -1, 0, x) - G{~l/y, 0, 0, x) 

+2G(-y,-l,0,a;)-G(-y,0,0,x) 

+2C2-ff (-1, x) + 4i/(0, a:) - 4C2i/(0, y) 

+G(-l/y, 0, x)iJ(0, y) + G(-y, 0, a;)i?(0, y) 

+2i/(0, a;)iJ(0, y) + 2G(-l/y, 0, x)H{l, y) 

+2G(-y, 0, x)H{\,y) + 4i/(0, a;)if (1, y) 

-4i/(-l, 0, a;) - 277(0, y)H{-l, 0, a:) 

-4i/(l, y)iJ(-l, 0, x) + 2i7(0, 0, x) 

+2H(0, a;)i?(0, 0, y) + 4iJ(0, a;)i?(0, 1, y) 

+4i/(0, x)H(l, 0, y) + 8if (0, x)H{l, 1, y) 

+4iJ(-l, -1, 0, x) - 2iJ(-l, 0, 0, x) 

-AH{0, -1, 0, x) + 2iJ(0, 0, 0, x) 



^H{0,0,0,y)-2H{0,0,l,y)-2C3 



+ Oie). (9) 



The 2-diniensional harmonic polylogarithms G 
are introduced in Appendix A of |24j. Here, we 
use the notations of |J31; the G(— y, 0,a::) and 
G{—y, — 1, 0, x) may be directly taken from (C.27) 
and (C.31). Further it is: 



G(-y,0,0,x) 



1 r dv In V 



2 Jo v + y 



In X 



ln(l + -) + lnxLi2(— ) - Li3( ). (10) 

2 y y y 

To summarize, the next, certainly ambitious 
task is the calculation of the remaining box type 
master integrals. Finally, for applications at a 
linear collider, the complete 2-loop QED am- 
plitudes, the 1-loop Standard Model corrections 
|25l26j ■ and cut dependent real photon radiation 
effects (determined with Monte Carlo generators) 
must be evaluated. 
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